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Abstract. Let V be a potential on R 3 that is smooth everywhere except at a discrete set 
S of points, where it has singularities of the form Z/p 2 , with p(x) = \x — p\ for x close to p 
and Z continuous on R 3 with Z(p) > —1/4 for p £ <S. Also assume that p and Z are smooth 
outside 5 and Z is smooth in polar coordinates around each singular point. We either assume 
that V is periodic or that the set S is finite and V extends to a smooth function on the radial 
compactification of R 3 that is bounded outside a compact set containing S. In the periodic 
case, we let A be the periodicity lattice and define T := R 3 /A. We obtain regularity results in 
weighted Sobolev space for the eigenfunctions of the Schrodinger- type operator H = — A + V 
acting on L 2 (T), as well as for the induced k-Hamiltonians obtained by restricting the 
action of H to Bloch waves. Under some additional assumptions, we extend these regularity 
and solvability results to the non-periodic case. We sketch some applications to approximation 
of eigenfunctions and eigenvalues that will be studied in more detail in a second paper. 
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1. Introduction and statement of main results 

We study in this paper regularity and decay properties of the eigenfunctions of Schrodinger 
type operators with inverse-square singularities. We either assume that the potential is periodic or 
that it has a nice behavior at infinity and only finitely many singularities. In order to explain our 
assumptions and results in more detail, we organize our Introduction in subsections, concentrating 
on the case of periodic potentials, the non-periodic case being similar, but simpler. We first 
introduce the operators obtained from the Hamiltonian —A + V acting on Bloch waves. In 
the second subsection of the Introduction, we explain our assumptions on the potential V. Finally, 
we state our main results and we summarize the contents of the paper. 

This paper is written to put on a solid foundations the numerical methods developed in [24] 
and [25] • We have thus written this paper with an eye to the numerical analyst. More theoretical 
results on Hamiltonians with inverse square potentials in arbitrary dimensions will be included 
in the third part of this paper. 

We have to mention Kato's ground breaking papers [30], where the self-adjointness of Schrodinger 
type Hamiltonians was proved and [3T], where boundedness properties of the eigenfunctions and 
eigenvalues of these Hamiltonian operators was proved. Moreover, see [SJ |H1 US HH EH HH1 EH] for 
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other papers studying Hamiltonians with inverse square potentials, both from the point of view 
of physical and numerical applications. See also [l[ini[Il[Il[I3[M[2ni[23lllSZl[Sn]for some 
related results. 

1.1. The Hamiltonian H^. Let V be a periodic potential on M 3 with Bravais lattice (of trans- 
lational symmetries) A = Z 3 . Assume that V is smooth except at a set of points S, which is 
thus necessarily also periodic with with respect to the lattice A. We assume that there are only 
finitely many elements of S in any fundamental domain P of A. Let p G S be a singular point 
and p(x) ~ \x — p\ for x close to p and p smooth outside iS. We assume that around p the 
potential V has a singularity of the form Zj p 2 , where Z is continuous across p and smooth in 
polar coordinates around p. We shall study numerically Hamiltonian operators of the form 

(1) H:=-A + V. 

Systems with such potentials have been studied as theoretical models both from the viewpoint 
of classical mechanics and from the quantum mechanical viewpoint. In addition, they arise in 
a variety of physical contexts, such as in relativistic quantum mechanics from the square of the 
Dirac operator coupled with an interaction potential, or from the interaction of a polar molecule 
with an electron [42] . 

A standard method for studying these operators is through their action on Bloch waves. For 
any kgl 3 . the Bloch waves of H with wave vector k are elements of L 2 oc (R 3 ) that satisfy the 
semi-periodicity condition that, for all I 6 A, 

(2) Mx + X) = e^ x ^{x). 

(It is enough to consider k in the first Brillouin zone P* of the reciprocal lattice to A. Also, the 
equality is that of two Lf oc functions, and hence it holds only almost everywhere in x.) A Bloch 
wave with wavevector k can be written as 

(3) Mx) = e* k - x u k (x) 

for a function u k that is truly periodic with respect to A t and thus can be considered as living 
on the three-torus T := K 3 /A ~ (S' 1 ) 3 (obtained by identifying points in K 3 that are equivalent 
under the action of the lattice A by translations). Note that the periodicity condition that a 
Bloch wave satisfies prevents it from being in L 2 (K 3 ), thus a nontrivial Bloch wave that satisfies 
the equation 

(4) tf^k = A^ k 

is not, properly speaking, an eigenfunction of the Hamiltonian operator H. Rather, it is a 
generalized eigenfunction, corresponding to a value in the continuous spectrum of H. If ipk is 
a Bloch wave that is a generalized eigenfunction of H with generalized eigenvalue A, then the 
function «k := e _lk ' :E ?/'k(a;) will then be an actual A-cigcnfunction of the k-Hamiltonian on 
L 2 (T) defined by 

3 

(5) H^:=-Y^(dj+ik.j) 2 + V. 
Indeed, this follows from the equation 

(6) H(e^ x u k (x)) = e^H^x). 

Thus, it is useful to understand the regularity of eigenfunctions Uk for the operators H^, as well 
as to arrive at theoretical estimates for the accuracy of various schemes to estimate them and 
their associated eigenvalues. 
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1.2. Assumptions on the potential V. In this paper, we extend and test the results of [37] 
to deal with the more singular potentials that have inverse-square singularities. More precisely, 
we extend the results of the aforementioned paper from potentials where pV is smooth in polar 
coordinates to potentials where p 2 V is smooth in polar coordinates and continuous on T. (Recall 
that p is a function that locally gives the distance to the singular point.) In particular, we obtain 
regularity results in weighted Sobolcv spaces that will then permit us to derive estimates for the 
accuracy of two approximation schemes that we design and which are studied in detail in the 
forthcoming second and fourth parts of this paper. The first scheme is a finite element method 
with a mesh graded towards the singular points as in [271 [7] . The second scheme is an augmented 
plane- wave method, similar to a "muffin-tin" method |38j . 

In order to state our results, we first need to set some notation and introduce our assumptions 
on the potential V. Let S C T := R 3 /A be the finite set of points where V has singularities. By 
abuse of notation, we shall denote by \x — y\ the induced distance between two points i,j £ T. 
Let then p : T — > [0, 1] be a nonnegative continuous function smooth outside S such that 

(7) p(x) = \x—p\ for a; close to p e S, 

as before, and further assume also that p{x) = 1 for x far from S. 

Our first assumption on V is that p 2 V be smooth in polar coordinates up to p = near each 
singularity. Let us explain this in more detail. We first replace each singular point p £ S with 
a 2-spherc in a smooth way, thus obtaining a space denoted T \ <S. This is the usual procedure 
of blowing up the singularities. We think of stretching out the holes where the singularities of 

V are and compactifying the result using boundary spheres. It would be possible to carry out 
analysis similar to the calculations in this paper with only the assumption that p 2 V be smooth on 
T \ iS. To simplify some calculations and to obtain closed form results, however, we will further 
require the resulting function Z to be constant on the blow up of each point in 5, which can be 
reformulated as saying that p 2 V is also continuous on T. Our first assumption on the potential 

V is therefore 

(8) Assumption 1 : Z := p 2 V G C°°(TT~S) n C(T). 

Assumption 1, more precisely the continuity of Z at S, allows us to formulate our second assump- 
tion. Namely, 

(9) Assumption 2 : rjo := minZ(p) > —1/4. 

pes 

Therefore, the constant 

(10) T) := y/T/I+m, 

which will play an important role in this paper, is a positive real number. This constant will 
appear in many results below. We will use Assumptions 1 and 2 throughout the paper, except in 
Section [2] where we prove more general forms of our results, not requiring Assumption 2. 

1.3. Regularity and approximation results. The domains of all the Hamiltonian operators 
considered in this paper will be contained in weighted Sobolev spaces on T \ 5. We define these 
spaces by: 

(11) /C™(T \ S) := {v : T \ S -> C, p^-'^v e L 2 (T), V < m}. 

These spaces have been considered in many other papers, most notably in Kondratiev's ground- 
breaking paper [32]. They can be identified with the b-Sobolev spaces of [39] (associated to a 
manifold with boundary), but with a different indexing and notation. These spaces were gener- 
alized in [3] to more general manifolds with corners with additional structure (Lie manifolds). 
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To formulate the stronger regularity for eigenvalues, we shall need the following notation. For 
each point p £ <S, let 

r 2 z{p) > | 

(12) MP) = { l + VW+Z(p) ^(P)G(-J'I) 
and 

(13) v Q = min^ (p). 

pes 

For each point p £ £ for which Z(p) £ ( — 1/4,3/4], define a smooth cutoff function Xp that is 
equal to 1 in a small neighborhood of p and is zero outside another small neighborhood of p. so 
that all the functions Xp have disjoint supports. Define the space W s to be the complex linear 
span: 

(14) W s = C XpP ^+^- 1/2 . 

Z(p)G(-l/4,3/4] 

Using also the notation introduced in the previous subsection, we then have the following result, 
whose proof follows from the proof of Theorem 12.31 below. 

Theorem 1.1. Consider a potential V satisfying Assumptions 1 and 2. Then the Hamiltonian 
operator acting as an unbounded operator on L 2 (T) has a distinguished self-adjoint extension 
with domain 

T>{H k )= K? 2 (T \ S) + W s C /C^(T \ S), v<v Q = min v (p) € (0, 2] . 

pes 

In particular, if r) := rjo := min p6l s Z(p) > 3/4, then is in fact essentially self-adjoint and, if 
f]a > 3/4, then D(if k ) = /C|(T \ «S). 

The importance of the above theorem is the following corollary, which says that under As- 
sumptions 1 and 2, the Hamiltonian operators = — Ylj=i(®j + ikj) 2 + V acting on L 2 (T) can 
be completely understood through their cigenfunctions and eigenvalues. 

Corollary 1.2. Under the assumptions of Theorem there exists a complete orthonormal 
basis of L 2 (T) consisting of eigenfunctions of . 

We can now state a regularity theorem for the eigenfunctions of near a point p £ S, or 
cquivalcntly, for Bloch waves associated to the wavevector k. Recall the functions Xp supported 
near points of S and used to define the spaces W s . 

Theorem 1.3. Assume that V satisfies Assumptions 1 and 2. Let H^u = Xu, where u £ P(iJk), 
Then, for any m £ Z+, 

u £ (T s S), Va < n := min y/l/4 + Z(p). 

pes 

Moreover, we can find constants a p £ R such that 

U -Y. X P P^ 1/4+Z{p) ~ 1/2 G IC n J+l (Ts5), Va' < min ^9/4 + Z{p) . 
pes peS 

The next result, which is the last we will mention in this introduction, will permit us to 
construct approximation schemes for the solutions of equations of the form (A + Hk)u = f. 

Theorem 1.4. Let us use the notation of Theorem \1.1\ and both Assumptions 1 and 2. Then 
there exists Co > such that A + f/k : ^CJJ+i^T \ S) — > K%~i (T \ S) is an isomorphism for all 
m £ Z>o, all \a\ < r], and all A > Cq. In particular, Hk is symmetric and bounded below, thus has 
a Friedrichs extension, which is equal to the closed extension considered in Theorem \l.l\ above. 
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From now on, we shall write iJk for the Fricdrichs extension of the original operator defined 
in Equation © and P(iJk) for its domain. 

We observe additionally that with the exception of Corollarv ll.2[ all of the results above extend 
to Hamiltonian operators on M. 3 associated to a non-periodic potential with a finite number of 
inverse square singularities satisfying Assumptions 1 and 2 and radial limits at infinity. This is 
because the techniques employed to obtain the results are local, and a Hamiltonian operator over 
R 3 with a smooth potential that has radial limits at infinity is always essentially self-adjoint. Of 
course, in this situation, there are only isolated eigenvalues below the continuous spectrum, and 
the bulk of the spectrum is continuous. 

The paper is organized as follows. In Section [5] we prove Theorem 11.11 identifying a closed self- 
adjoint extension of the operator H^, and Theorem 11.31 giving regularity results for cigenfunctions 
of this closure. Some of the results of this section do not rely on Assumption 2. Beginning with 
Section |3l however, we shall require that Assumption 2 be satisfied. In that section, we prove 
that for rj > —1/4, Hk is bounded from below, and we can thus identify the closure from Section 
[2] as the Friedrichs extension of H^. In Section 0] we discuss how our results extend to the the 
nonperiodic case and how to use them in numerical methods. 

Acknowledgements. We would like to thank Alexander Strohmaier, Joerg Seiler, Thormas 
Kraincr, Jorge Sofo, and Anna Mazzucato for useful discussions. We also thank the Levcrhulmc 
Trust whose funding supported the fourth author during this project. This project was started 
while Hunsicker and Nistor were visiting the Max Planck Institute for Mathematics in Bonn, 
Germany, and we are grateful for its support. 

2. Regularity and singular values 

The regularity analysis of the operators is done locally in the neighborhood of each pe5. 
Let us recall that Z := p 2 V G C°°(T \ <S)nC(T), which is our Assumption 1, which we will require 
to hold true throughout this paper. In this section, we shall mention explicitly when Assumption 
2 is used, since some of the results hold in greater generality. 

For simplicity of the notation, we shall assume that S consists of a single point p. The results for 
potentials with several singularities with different values of Z(p) can then be pieced together from 
local versions of the result in the one singularity case. The proofs of Theorems 11.11 and Theorem 
11.31 (which we prove in this section in more general forms not requiring Assumption 2), rely on 
the pseudodifferential operator techniques of the b-calculus and b-operators [4j [21] [39l [36l 05] . A 
review of these basic tools is contained in [27j , so we will not go into detail about them again here. 
Throughout this paper, we will refer to b-operators and the b-calculus, although the properties 
can be equivalcntly described in terms of cone operators and the cone calculus, and in fact, in 
some of the references in this section, they are referred to in this way. For a discussion of the 
equivalence of the b- and cone calculi, see [35] . 

2.1. The boundary spectral set. In order to use the b-calculus, we study the associated b- 
diffcrcntial operators 

We can write such an operator in polar coordinates around p G S as 
(15) P k ,A = (pd P f + pd p + As»-i - p 2 V - P B KX , 

where 

B K x :=P ($^(-2*0, + ^)^ 

is a first order b-operator. 

The operator Pk,A is an elliptic b-opcrator on T \ <S. We calculate the indicial family of Pk,A 
at a point p G S, denoted (Pk.\) p (T), by replacing pd p with r in Equation (|15j) . and by replacing 
the coefficients with their values at p = 0. Doing this, we find that the indicial families for Pk,A 
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at p G S are, in fact, independent of A and k, since the dependence on k and A affects only B^ x, 
and p vanishes at the singular points. We summarize this discussion in the following lemma. 

Lemma 2.1. Let P := Po,o- Then the indicial family of Pk.A at p € S is given by 

(16) (P^x) P (r) = P p (r) = r 2 + t + A s2 - Z(p). 

Lemma [2.11 allows us to calculate the boundary spectral set Spcc b (P p ) for Pk,A at a given p. 
The boundary spectral set for Pk is then defined by ( |39j ) 

Spec b (P p ) := {(t, n) | P p (t) _1 has a pole of order n+ 1 at r}. 

By Lemma |2~T1 the set Spec b (P p ) will also be independent of k and A. To calculate the spectral 
set Spec b (P p ) explicitly, recall that the eigenvalues of A52 are — 1(1 + 1), for I G Z>q, and define 



h7 s R . y/(l + 2/)2+4Z(ri-l 

(17) Pl, P ■= ^ , / G 



(18) a;, p := , I G Z; 



By an abuse of notation, we take -^/(l + 2/) 2 + AZ(p) to denote the positive imaginary root when 
the quantity under the root is negative. Our discussion gives the following. 

Lemma 2.2. If Z(p) £ {-(1/2 + l) 2 }^, we have 

(19) Spec b (P p )= (J {GSi.p.Oj.Kp.O)}, 

;gz> 

anrf, «/ ^(p) = —(1/2 + l p ) 2 , for some l p > 0, l p G Z, 

(20) ^ec b (P p ) = {(-l/2,l)}U |J {(A, p ,0),K Pl 0)}. 

ZGZ> 

When Z(p) < —1/4, a finite number of these values of a and /3 will be complex with real part 
equal to —1/2. This is one of the reasons why we have introduced Assumption 2, which states 
that Z(p) > —1/4 for all p G S. Of course, if Assumption 2 is satisfied, Spec b (P p ) is given only 
by Equation (fT9|) . The case when Z(p) is close to —1/4 is important because it gives to some 
interesting numerical phenomena and also because in various applications, it has an interesting 
interpretation, see for instance [4"2"] . 

The machinery of the b-calculus now gives us information about closed self-adjoint extensions 
of .Hk (see [36l [19\ HH 04] for details). Note that in this setting, we are considering extensions of 
our operator acting on the core consisting of smooth functions supported away from the points 
of S. For Z(p) < 3/4, there will be several possible self-adjoint extensions. Compare this to the 
case of extending the Laplacian operator on T from acting on the core consisting of all smooth 
functions. In this case, there is a unique self-adjoint extension. This is because the core is larger 
than in our case. The extension obtained using the larger core is one of the possible extensions 
obtained using the smaller core, but it is not the only one. This is why one does not see the issue 
of choosing a sclf-adjoint extension arising when the potential is of the form Zp a for a > —2, for 
instance, in the Coulomb case considered in |27j . 

With the above lemmas in place, we can now prove Theorem 11.11 In fact, we shall prove a 
stronger result that does not require Assumption 2. 

Theorem 2.3. Consider a potential V satisfying Assumption 1 and assume that S consists of 
just one point p. Then the Hamiltonian operator acting as an unbounded operator on L 2 (T) 
has distinguished self-adjoint extension with domain 2?(iJk) C IC 2 (T \ S) for all v < G (0,2]. 
In particular, if Z(p) > 3/4, then is in fact essentially self-adjoint and, if Z(p) > 3/4, then 
"D(Hk) = A^2 (TT \ S). If Assumption 2 is satisfied, we also have 



V{H^)=K 2 2 {T\S)+<C X p r '- l/ \ r) = VW+Z(pj- 
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where x is a cutoff function that is zero outside some neighborhood of p and equals 1 close to p. 

Proof. For each k and A, the operator — A is a symmetric, unbounded b-operator on L 2 (T) 
(see [39l EH US]). Define the operator A = p 1 / 2 H\ L p~ 1 / 2 . Then A is a symmetric unbounded 
b-operator on p~ 1 Ll(T \ S) — /C°y 2 (T s S). The self-adjoint extensions of A correspond exactly 

to those of -f/k — A with domains shifted by weight p 1 ^ 2 , so we will study the self-adjoint extensions 
of A as the calculations are somewhat easier in this case. 

By Lemma [2TTT the indicial roots of A are the roots of shifted by 1/2. So we let 



(21) fa = y/(l + 1/2) 2 + Z{p) and &i = -fa. 

Note that 0/ftel HZ(p) > -(Z + ±) 2 , there is a double root at fa = if Z{jp) = -(Z + ±) 2 , and 
^ fa £ iM. if Z(p) < — (I + \) 2 ■ The critical strip for self-adjointness of unbounded operators 
on p~ l L\{^ s S) is (—1,1), that is, an operator is essentially self-adjoint if and only if it has 
no indicial roots with real part in this interval (see (36] [19]). Recalling that I £ Z>o, we see 
that for Zip) > §, there are no roots in the critical strip so the operator is essentially self- 
adjoint. If further Z(p) > | we get the somewhat stronger result that 2?(iJk) = /C 2 (T \ S). For 
Z(p) £ (—\, |) we get two real roots in the critical strip corresponding to I = 0, for Z(p) = — j, 
we get a double root at in the critical strip corresponding to I = 0, and for Z{p) < —j, we get a 
finite number of complex conjugate imaginary root pairs and possibly two real roots or a double 
root at in the critical strip corresponding to some finite set of I. 

By the theory in |19j . the space £ := Dom(A max )/Dom(A m i n ) is finite dimensional and spanned 
by functions local around p of the form: 

i I 

(22) u u {V'ff,«)p-Aff}u u u {<,«ih^n, 

|SR(A)|e(o,i)™=-i ft=o™=-< 

where w; is a local cutoff function that equals 1 near p and for p large, and where the ip™ are an 
orthonormal basis for spherical harmonics with eigenvalue 1(1 + 1). Further, the operator A with 
domain T> := Dom(A m ; n ) + span(ui, . . . , u n ) is self adjoint if, and only if, linear combinations of 
these basis functions form a maximal set on which the pairing, [it, vJa is trivial, where 

(23) [u, v]a ■= j> Au(a) - s i v(a) da. 

Here 7 is a simple closed loop around the indicial roots of A in the critical strip, " represents the 
Mcllin transform and -52 denotes the standard L 2 paring on S 2 . Since the tp™ arc orthonormal, 
this pairing reduces to a sum of loop integrals of the form: 

(24) [u h Vi]A = ~^f(° 2 +ff)ui(o-)i^)do- 1 

where ui = u+wp^ 1 + u-wp~P' and vi = v+wp^ 1 + v-wp~@ l if fa ^ and ui = u+w + U-tvlog(p) 
and vi = v + w + v-wlog(p) if fa = 0. 

We can consider three cases: fa > 0, fa = and fa £ iR. First, as in [E5], define 

/>oo 

= -pd p w(a) := - / p-'° w' \p)dp. 
Jo 



Then we get: ^>(cr) = $(— a) and $(0) = 1. Also, using the properties of the Mellin transform, 
we find that for any a £ C, 

— TT. $(<r ± ifa) 
wp±Pi(<r) = — 5— ■ 

Now consider the case fa > 0. Carrying out the loop integral by evaluating residues, we arrive 
at the equation 

[ui,vi] A = fc(lt + U_ - lt_W+) 



8 



E. HUNSICKER, H. LI, V. NISTOR, AND V. USKI 



for a constant k ^ 0. If we set = 0, this reduces to arg(u + ) = arg(u_). Thus to get a 

sclf-adjoint boundary condition at p for A we can fix any ratio of \u + \ to We will choose to 
enlarge the minimal domain by the set with \v,-\ = 0, so spanned by {wp^ l ipj n } l m= _ l . Note that 
if Zip) > —j, we have / = 0, so to create a self-adjoint extension of A = p 1 ! 2 H-^p^ 1 / 2 , we only 

need to expand the minimal domain by the span of wp q = w/)\/ 1 / 4+z W. 
Next consider the case when (3i = 0. In this case we get 

- §((7) *'(<r) 

w{<7) = and wiogp = — » • 

a <j z a 

Carrying out the loop integral in this case, we arrive at the equation 

[ui, vi]a = fc(u + v„ + U-V + ). 

Again setting [ui, ui]a = 0, we arrive this time at the condition it + u_ € iR. So we may this time 
again choose to fix w_ = and enlarge the minimal domain by the set spanned by {wipl a } l m= _ l . 
Finally, consider the case when Pi = ia. This time we get the equation 

[ui,vr) A = k(u + v + 

Setting [ui,u{\ = we arrive at the condition = We can choose u + = 1 and 

u_ = — 1 to get a self-adjoint condition by enlarging the minimal domain by the set spanned 
by {wcos(alogp)^ ; m }J n= _ r 

In order to get back to the corresponding choice of self-adjoint extension of — A, we multiply 
each basis element by p~ x l 2 . Since each of these basis functions is in p M L 2 (T) for all fj, < 1, overall 
we find that V{H) C p^L 2 (T) for all \x < 1. This completes the proof of ThcoremO □ 



2.2. Some corollaries. We now prove some consequences of Theorem l2.3l First, its proof implies 
the following stronger result: 

Corollary 2.4. If r] > 0, then V(H k ) C H 2 (T) n p e C°(T), for some e > 0. 



We now deduce Corollary 1 1 . 2 1 from Theorem 11.11 and its proof. This corollary is specific to the 
periodic case. 

Proof, (of Corollary 1 1.2[) . The extension of is self-adjoint, hence has real spectrum, so we get 
that for A ^ R, the operator 

A - H k : V(H k ) -> L 2 (T) 

is a bounded invertible operator with bounded inverse Qk,A, which is the resolvent of at A. By 
Theorem 11.11 by the definition of the weighted Sobolcv spaces /C™ (T \ S) , and by the b-Rcllich 
lemma [3 [39] , we have for some e > that 

V(H k ) c /C 2 (T\5) <e£ 2 (T). 

(Recall that (s means "compactly embedded".) Thus for A ^ R, the resolvent of i?k, 

ftx(ffk) : L 2 (T) -> L 2 (T) 

is a compact operator. By standard results of functional analysis, if a self-adjoint operator has 
compact resolvent, then L 2 has a complete orthonormal basis consisting of eigenfunctions for this 
operator. □ 

2.3. Singular functions expansion. To prove Theorem 11.31 we again use results of the b- 
calculus, this time primarily from |39j . Again, we prove a more general statement that does not 
require the Assumption 2. Let vq be as in Equation (|12l) . 



Theorem 2.5. Assume S = {p} and Z := p 2 V satisfies Assumption 1. Assume H^u = Xu where 
u £ 2?(iJk)- Then for any m G Z + and any v < vq, 

ue/C™(T\<S). 
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Further, near each p G S, where Z(p) ^ (I + 1/2) 2 for any I > 0, u has a complete (though not 
unique) expansion of the form: 

(25) u = u 

where the formula for 1z(p) is given in Equation (|27jl below, uo is smooth up to p — in polar 
coordinates and vanishes to all orders there, hence is in fact smooth on T and vanishes to all 
orders at p, and the coefficient functions <? 7 are smooth functions on S 2 . Under the additional 
Assumption 2, when Zip) > —1/4, the first coefficient g v —i/2 is a constant function. 

Proof. Any cigcnfunction u of must be in its domain, thus in 1C 2 (T \ S) for all v < vq. Our 
first goal is to improve the degree of smoothness from 2 to m for m G N. To do this, we use the 
fact that any A G C, the operator — A is Fredholm as a map between weighted Sobolev spaces: 

(26) ff k - A : K£ (T \ S) -> JC n a l S 2 (T x S) 

for all a G K such that a £ U; e z >0 {/3; iP + i,cti,p + §} = Spec b (i?k) + 3/2. By general b-calculus 
theory, the set Spec fc (i/k) is a discrete subset of C and furthermore, for any 70 and 77, it has only 
a finite number of elements in the strip 70 < $t(z) < 77. Thus for any uq, there exist arbitrarily 
close v < Vq such that the condition on a is satisfied for a = v + 2s, where s G N. Together with 
standard bootstrapping arguments, this allow us to improve the regularity of eigenfunctions of 
fljk in terms of weighted Sobolev spaces to IC" l (T \ S) for all m and v < i/q. 

Next, to obtain the expansion in 11.31 we use a general result in the b-calculus literature, see 
e.g. [39) . that implies that any u G L) m . a IC" l (T % S) which is an cigcnfunction for in some 
weighted L 2 (T) in fact has much stronger regularity: it polyhomogeneous in p near each p G S 
with index set Iz( P )- If Z(p) ^ { — (I + ^) 2 }^Z (for instance if Z{p) > —1/4), then the index set 
is simply a set of complex numbers that is finite in any strip 70 < 5ft(z) < n: 

00 

(27) Iz( P ) = U {Pl, P + n, a liP + n} l£Z > . 

n=0 

This means that around each p G S, there exist smooth coefficient functions g 1 G C°°(S 2 ) such 
that for all N, 

(28) ^^^-^^ep^TKS), 7 € l Z ( P )^(l) < N 

that is, un is iV times continuously diffcrcntiablc up to p = in polar coordinates at p and in 
these coordinates, vanishes with all of its derivatives up to order p N there. We take the limit in 
the topology of the smallest /C™ to which u belongs. 

If Z(j>) ^ -(I + \) 2 for any / > and if u G L 2 (T), when we let N -> 00, we find that 

« = + 5> 7 <? 7 , 7€I,5R( 7 ) > -1/2, 

where uq is smooth up to p — in polar coordinates and vanishes to all orders there, hence is in 
fact smooth on T and vanishes to all orders at p. 

Since the set of 7 that appear in this expansion is discrete in K, we get that the smallest 
exponent that appears will in fact be somewhat better than —3/2. This first exponent will be 
if Z(p) < 3/4. If Zip) > 3/4, then eigenfunctions will in fact be in a space with higher weight 
than /^(T \ S): the weight will be 1 + Zip) + 1/4. 

Finally, we can note that the terms of the expansion of an eigenfunction for that are not in 
/C|(T \ S) will be of the forms determined in the proof of Thcorem ll.il (see, eg [55] for a proof). 
So, for instance, if Zip) > —1/4, the leading term of any eigenfunction u will be constant in S 2 , 
and u minus its leading term will vanish at p. Further, if Zip) > —1/4, then the exponents 7 will 
all be real numbers. Thus we obtain: 

(29) u = uo + p n ^^g r ,-i/2 + ^P 1 9 1 , 7 e x z( P )l > V - V 2 , 

where g^-i/i a constant. This completes the proof of Theorem II .31 □ 



10 



E. HUNSICKER, H. LI, V. NISTOR, AND V. USKI 



3. Invertibility 

In this section we prove the boundedness and invertibility result in Theorem 11.41 From now 
on, we require both Assumptions 1 and 2 to be satisfied by our potential V . 

3.1. Preliminary results. We begin with a few standard results lemma. 

Lemma 3.1. Let ro,ael. Then 
(i) For any f G C°°(T \ S), the multiplication map 

£™(T \5)3ii4/u€ /C™(T s S) = p a - 3 / 2 H^(T \ S) 

is continuous for all m £ Z + and all a 6l. 
(ii) The operator — A maps /C^S (T \ S) to fc™Li (T \ S) continuously. 
(Hi) The operator p B^x maps /C^+^T \ S) to /C™(T \ S) continuously, 
(iv) p- l B KX : JC+i 1 ^ \ ^) "> ^a-Y(T \ 5) is compact 

Proof. The simple proofs of these results are the same as that of the analogous results in [27] . 
and follow directly from properties of b-operators [3] [39l [36] . □ 

We also need the following standard lemma (again, see [STJ for its proof). 

Lemma 3.2. Let aelle arbitrary and assume that u £ /Cf +a (T \ 5) and that v £ /C^_ a (T\<S). 
Tfterc (Aw, v) + (Vu, Vu) = 0. 

We shall also need the following consequence of the general properties of the b-calculus [39] 05] . 

Proposition 3.3. Let us fix A £ C and a ^ {A, p ,az, p } = U; e z >0 {/3/ iP + |,a; lP + |}. -Let 6e 
£/ie number of elements in the set {/3/,p, ck;, p } i/iai are between ana 1 a, counted with multiplicity. 
Then the operator i/k — A is Fredholm as a map between weighted Sobolev spaces: 

ff k - A : /C™V (T \ <S) -> /CM 1 (T x S) 

and has index —N if a > 0, respectively N if a < 0. 

Proof. We consider again the operator Po,o = p(-ffk — A)p, which is a b-differential operator. It 
is unitarily equivalent to p^-Po.oP -1 ^ 2 acting on b-Sobolev spaces (see the proof of Theorem 
12. 3p . which has {/3z,p, a;.p} as a b-spectrum. The result then follows from the characterization of 
Fredholm b-differential operators [511 1331 EH] • 

It remains to determine the index of — A. Let m a be the index of the operator for a fixed 
value of a. Then it is a standard result that m a — mb is given by the number of singular functions 
with exponent between a and b [33] [39] 05] [41] . This is enough to complete the proof. □ 

See [25] for an extension of this result and for more details. 
Now recall the Hardy inequality, which states that 

(30) f r- 2 \u{x)\ 2 dx< {2/{N-2)) 2 f \Vu(x)\ 2 dx 

Jm N Jm n 

for any u £ H 1 (R N ), N > 3, where r is the distance to the origin [17]. We can use this to prove 
the following important lemma. To simplify notation, after the lemma statement, we shall let 
(u,v) := (u,v) L 2 {T) . 

Lemma 3.4. There are constants C, 7 > such that for any u £ K.\(T \ S), 
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Proof. For an operator T : K.\(T \ S) -» /Cl}(T \ 5), we shall write T > if (Tu,u) > 
for all u G K.\(T \ <S). Now let > be a smooth function T that is equal to 1 in a small 
neighborhood of S and has support on the set where p(x) is given by the distance to S and let 
Vo(x) = Z{p)4>(x)p~ 2 (x) for x in the support of 4> and close to p G S. Outside the support of </>, 
we let Vq = 0. Then Hardy's inequality applied to (j^^u, which we can think of as living now on 
M 3 rather than T, gives 

(31) (<j) 1/2 {-A + zV )<j) 1/2 u,u) > and (^ 1 / 2 (-A)0 1 / 2 w, u) > 0. 

We can think of this as saying that the most singular part of the operator H k , that is, T = 
0i/2 (_ A + zV )(f> 1/2 , satisfies T > 0. We will prove Lemma 13.41 by decomposing the operator 
i?k + C as a sum of four operators 

-^k = Tl + T 2 + I3.C + ^4,(7) 

which we will show are all bounded from below for sufficiently large C. 
Recall we can write 

H k = -A + V + V 1 +p- 1 B k .o, 

where V± := V - V satisfies pVi G C(T) and p -1 ^. is a first order differential operator over T 
with smooth coefficients. 

Assumption 2 and Equation (|31j) imply that for e < 1, the operator 

(32) T x := (1 - e)T > 0. 

Fix any suitable value for e > 0. Then we can write + C in terms of T\ by decomposing in 
terms of the multiplication operators 1 / 2 and (1 — 1//2 : 

(33) H k + C = eT + T x - ^ 1/2 A^/ 2 - (1 - 1 / 2 )A(1 - 4> 1/2 ) +V X + R x , 

where R\ is a first order differential operator with smooth coefficients and ip = 2(f) 1 / 2 (l — (p 1 ^ 2 )- 
Let T 2 := -ij/^Aip 1 / 2 - (1 - 1 / 2 )A(1 - <f>^ 2 ). Then T 2 > by the Hardy equality applied to 
V> 1/2 w and to (1 - (^l 2 )u. Define T 3 := -eA + i?i + C/2 and T 4 = e<Pp~ 2 + V\ + C/2. We claim 
that for C large enough, T 3 > and T 4 > 0, which will prove the result. 

The proof that T 4 > for C >> follows from a straightforward calculation minimizing the 
function e<frp~ 2 + V\. The proof that T3 > for C » is basically the same as the proof that 
a Schrodinger operator with periodic Coulumb type potential is bounded below. This is proved, 
for example, in [27]. □ 

Note that the above lemma implies that H k is bounded from below as an operator /C}(T\5) — > 
/CI 1 (T \ S), which is the special case of Theorem 11.41 when m = a = 0. In addition, if we define 
the form a(u,v) := ((H k + C)u,v), where the right-hand side is the natural pairing between 
elements of JCz\ and K,\, then this lemma implies that a(u, v) satisfies the assumptions of the 
Lax-Milgram lemma for the vector space V = K,{(T \ S). This and Cea's lemma imply that for 
any element u £ JC\(T \ S) and any finite dimensional subspace V C /Q(T \ S) we can construct 
a unique (Galerkin) approximation uy S V for u that, up to a multiple independent of it, is the 
best approximation for u in V. 

If we could also use the K\ norm in our approximation results, we would now have the necessary 
tools to prove it. However, we need to use the slightly smaller space K} a+1 instead. Thus we the 
stronger result, Theorem 11.41 to ensure the Lax-Milgram theorem and Cea's lemma apply to the 
analogous form on these spaces. 

We shall also need the following regularity result. 

Proposition 3.5. Let a, A € M., m G Z+. There exists a constant C > such if u G K,\ +1 (T \ <S) 
and (A + H k )u G /C^T^T \ S) then u G ^^^(T \ S) and 

IMIjc-Y <C(\\{\ + H k )u\\ K ^ + \\u\\ KUi ). 
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Proof. We consider again the operator Po.o = p(Hk ~ A)p, which is a b-differcntial operator. Our 
result then follows from the regularity for b-pseudodifferential operators [35] 13] • O 

We now complete the proof of Theorem 11.41 as follows. 

Proof, (of Theorem ll.4[) . As in [37], by regularity for b-diffcrential operators, if we prove our 
result for m = 0, then the regularity result of Proposition 13.51 implies it for all m > 0. We shall 
thus assume m = and focus on extending Lemma |3.4[ where a = 0, to the case when \a\ < r/. 

Fix C as in LemmaGOO Let D a := C + i? k : ££ +1 (T \ <S) -> /C^ X (T \ 5), that is to say, 
C + -ffk with fixed domain and range. As usual, we may identify the dual of K\ (T \ S) with 
the space /C~^(T \ S) using the £ 2 -inner product. Then using the symmetry of H^, we find that 
D* a = D_ a . " 

By Lemma T3.41 the operator Do is invertiblc. By basic results in b-calculus, D a is Fredholm 
for \a\ < rj since the weighted spaces in its domain and range do not correspond to an indicial 
root as calculated in the previous section (see Proposition I3.3|) . Hence, for such a, the family 
p a D a p~ a is a continuous family of Fredholm operators between the same pair of spaces. Since 
index is constant over such families, we have that md(D a ) = for all < a < rj. We want to 
know these operators are all isomorphisms. By the index calculation, it now suffices to show they 
are all injective. 

The inclusion K,^ +1 (T \ S) C K\ (T \ S) allows us to compute (D a u, u) = (Vm, Vu) + (u, u) for 
u € /C^ +1 (T \ S), by Lemma \'S. 2 1 Assume D a u — 0, then (D a u,u) — 0, and hence u — 0. This 
implies that the operator D a is injective for < a < rj. Since it is Fredholm of index zero, it is 
also an isomorphism. This proves our result for < a < rj. To prove the result for — ?/ < a < 0, 
we take adjoints and use D a = (D- a )* . 

By the characterization in |19j of the Friedrichs extension of a b-opcrator which is bounded 
below, we can see that the extension we constructed in Theorem 11.11 is in fact the Friedrichs 
extension of H^. The proof of Theorem 1 1.41 is now complete. □ 

The fact that the domain of the Friedrichs extension is (C — H-k)^ 1 (L 2 (T \ S)) and the theorem 
we have just proved give us a second way to identify the domain of the Friedrichs extension of 
H^. Following the method of [27J, we see that when is Fredholm on /C|(T \ S), the domain 
of the Friedrichs extension of consists of the span of KL\ (T \ S) and of the singular functions 
that are in IC{(T \ S) but are not in /C^(T \ S). This can be used to obtain an alternative proof 
of Theorem 1 1.1 1 if V satisfies both Assumptions 1 and 2, as follows. 

Proposition 3.6. Let Cq be as in Theorem \1.4\ and W s be as in Equation ()14j) . Assume the set 
{/3/ iP ,a; )P } does not contain 1. Then for A > Cq, the map 

\ + H k : /C" +1 (T n S) + W s -> iC^- x (T \ S) 

is an isomorphism. 

Proof. Let T := A + -f/k with the indicated domain and codomain. Proposition 13. 31 shows that T 
is Fredholm with index zero. Since 

/C™ +1 (T \ S) + W s C /Ci(T s S), 

Theorem 1 1 . 41 shows that T is injective. Hence it is also surjective, hence an isomorphism. □ 

4. Extensions and numerical tests 

We now discuss the extension to the non-compact case and indicate some applications to 
numerical methods. 
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4.1. The non-compact case. Most of our results in the previous sections extend to the non- 
compact case. Let R 3 ra d be the radial compactification of M 3 . We assume that the set of singular 
points iS C R 3 is finite and we replace each of the points in S with a two-sphere (that is, we blow 
up the singular points). Let (K 3 \ S) rad denote the resulting compact manifold with boundary. 
By p we denote a continuous function p : M 3 — > [0, 1] that is smooth outside S, close to each p E S 
it has the form p(x) = \x — p\, and it is constant equal to 1 outside a compact set. (Thus the 
difference with the function p considered in the periodic case is that now p is constant equal to 
1 in a neighborhood of infinity.) Then in the non-compact case, our Assumption I on Z := p 2 V 
is replaced with 

(34) Assumption l' : Z := p 2 V e C°°((R 3 \ S) rad ) n C(K 3 ). 

Assumption 2 remains unchanged. 

We consider now H = — A + V instead of the restrictions H^. Assumptions 1' and 2 allow us to 
extend to H all the results for of the previous sections, except Corollary 11.21 and Proposition 
13.31 The weighted Sobolev spaces /C"'((]R 3 \ S) rad ) are defined in the same way (but using the 
new function p). 

Let b e the infimum of V on the sphere at infinity. Then Corollary 11.21 must be replaced with 
the following characterization of the essential spectrum a e (H) of H: 

(35) a e (H) = [b e ,oo). 

To prove this result, one needs also the Fredholm conditions for operators in the scattering or SG 
calculus [34] [40] [44] [43] . Then in Proposition ^. 31 one has to take A < b e . Of course, in Theorem 
11.41 one will have Co > — b e . 

However, in the non-compact case, for applications to numerical methods, our results on eigen- 
values and eigenfunctions must be complemented by decay properties at infinity. The following 
is proved as in [2], Theorem 4.4. See also JT} [28l [29] . Let r : R 3 — > M. be a smooth function such 
r(x) = |.t| for x outside a compact set. 

Theorem 4.1. Let V be a potential satisfying Assumptions 1' and 2. Also, let0<e< V(x) — A 
for x outside a compact set and be u an eigenvector of corresponding to A. Then e er u € 

Under Assumptions 1' and 2, a perturbation argument further yields following result on the 
decay properties of the eigenfunctions and the solutions of the equation (C + H^u = f. 

Theorem 4.2. Let us assume that C + H k : K,\{(R 3 \ S) rad ) /Cli((IR 3 \ S) rad ) is invertible 
(which is the case if C > Co, with Co as in Theorem ] 1.4}) , then for \a\ and |e| small 

is again invertible. 

Proof. The proof uses the same arguments used in the proof of Theorem II. 4[ the continuity of 
the family e er i/ke" er in e, and the regularity result l3~5l □ 

4.2. Applications and numerical tests. Let u be an eigenvector of i/k or the solutions of 
equations of the form (A + H^)u = f , with / smooth enough. Our results give smoothness 
properties for u. They also give decay properties of u in the non-periodic case. These properties, 
in turn, can be used to obtain approximation properties of u. Standard numerical methods results 
(Cea's lemma or the results reviewed in [5]) then lead to error estimates in the Finite Element 
Method for the numerical solutions of the equation (C + H^u = f or for the eigenfunctions of 
i?k- We have tested these approximation results in the periodic case using, first, a graded mesh 
and, second, augmented plane waves. In both cases, the tests are in good agreement with our 
theoretical results. These numerical and the needed approximation results will be discussed in 
full detail in the second and fourth parts of our paper [MJ [55] . 



14 e. hunsicker, h. li, v. nistor, and v. uski 

References 

[1] S. Agmon. Lectures on exponential decay of solutions of second- order elliptic equations: bounds on eigen- 
functions of N-body Schrddinger operators, volume 29 of Mathematical Notes. Princeton University Press, 
Princeton, NJ, 1982. 

[2] B. Ammann, C. Carvalho, and V. Nistor. Regularity for eigenfunctions of Schrodinger operators. Lett. Math. 
Phys. 2012. 

[3] B. Ammann, A. Ionescu, and V. Nistor. Sobolev spaces on Lie manifolds and regularity for polyhedral domains. 

Documenta Math (electronic), 11:161-206, 2006. 
[4] B. Ammann, R. Lauter, and V. Nistor. Pseudodiffcrcntial operators on manifolds with a Lie structure at 

infinity. Ann. of Math. (2), 165(3):717-747, 2007. 
[5] D. Arroyo, A. Bcspalov, and N. Heuer. On the finite element method for elliptic problems with degenerate 

and singular coefficients. Math. Comp., 76(258):509-537 (electronic), 2007. 
[6] I. Babuska and J. Osborn. Eigenvalue problems. In Handbook of numerical analysis, Vol. II, Handb. Numcr. 

Anal., II, pages 641-787. North-Holland, Amsterdam, 1991. 
[7] C. Bacuta, V. Nistor, and L. Zikatanov. Improving the rate of convergence of high-order finite elements on 

polyhcdra. II. Mesh refinements and interpolation. Numer. Fund. Anal. Optim., 28(7-8):775-824, 2007. 
[8] S. Bidwell, M. Hassell, and C.R. Wcstphal. A weighted least squares finite clement method for elliptic problems 

with degenerate and singular coefficients. Math. Comp. (in press). 
[9] M. Dauge. Elliptic boundary value problems on corner domains, volume 1341 of Lecture Notes in Mathematics. 
Springer- Vcrlag, Berlin, 1988. Smoothness and asymptotics of solutions. 
[10] V. Felli, A. Fcrrero, and S. Terracini. Asymptotic behavior of solutions to Schrodinger equations near an 

isolated singularity of the electromagnetic potential, preprint. 
[11] V. Felli, A. Fcrrero, and S. Terracini. Asymptotic behavior of solutions to Schrodinger equations near an 

isolated singularity of the electromagnetic potential. J. Eur. Math. Soc. (JEMS), 13(1):119-174, 2011. 
[12] V. Felli, E. Marchini, and S. Terracini. On the behavior of solutions to Schrodinger equations with dipole type 

potentials near the singularity. Discrete Contin. Dyn. Syst., 21(1):91— 119, 2008. 
[13] H.-J. Flad, W. Hackbusch, and R„ Schneider. Best Af-term approximation in electronic structure calculations. 

II. Jastrow factors. M2AN Math. Model. Numer. Anal, 41(2):261-279, 2007. 
[14] H.-J. Flad and Harutyunyan. Ellipticity of quantum mechanical hamiltonians in the edge algebra, preprint 
larXiv:1103.0207H . 

[15] H.-J. Flad, R. Schneider, and B.-W. Schulze. Asymptotic regularity of solutions to Hartree-Fock equations 

with Coulomb potential. Math. Methods Appl. Sci., 31(18):2172-2201, 2008. 
[16] S. Fournais, M. Hoffmann-Ostcnhof, T. Hoffmann-Ostenhof, and T. Ostergaard S0rensen. Analytic structure 

of solutions to multiconfiguration equations. J. Phys. A, 42(31):315208, 11, 2009. 
[17] F. Gazzola, H.-C. Grunau, and E. Mitidieri. Hardy inequalities with optimal constants and remainder terms. 

Trans. Amer. Math. Soc, 356(6):2149-2168 (electronic), 2004. 
[18] J. Gil, T. Krainer, and G. Mendoza. Geometry and spectra of closed extensions of elliptic cone operators. 

Canad. J. Math., 59(4):742-794, 2007. 
[19] J. Gil and G. Mendoza. Adjoints of elliptic cone operators. Amer. J. Math., 125(2):357-408, 2003. 
[20] M. Griebel and J. Hamackers. Tensor product multiscale many-particle spaces with finite-order weights for 

the electronic Schrodinger equation. Zeitschrift fur Physikalische Chemie, 224:527-543, 2010. 
[21] D. Grieser. Basics of the 6-calculus. In Approaches to singular analysis (Berlin, 1999), volume 125 of Oper. 

Theory Adv. Appl, pages 30-84. Birkhauser, Basel, 2001. 
[22] T. Hausel, E. Hunsicker, and R. Mazzeo. Hodge cohomology of gravitational instantons. Duke Math. J., 

122(3):485-548, 2004. 

[23] B. Helffer and H. Siedentop. Regularization of atomic Schrodinger operators with magnetic field. Math. Z., 
218(3):427-437, 1995. 

[24] E. Hunsicker, H. Li, V. Nistor, and V. Uski. Analysis of Schrodinger operators with inverse square potentials 

II: FEM and approximation of eigenfunctions in the periodic case, preprint, submitted for publication. 
[25] E. Hunsicker, H. Li, V. Nistor, and V. Uski. Analysis of Schrodinger operators with inverse square potentials 

III: the n dimensional case. Work in progress, tentative title. 
[26] E. Hunsicker, H. Li, V. Nistor, and V. Uski. Analysis of Schrodinger operators with inverse square potentials 

IV: plane waves and orbital functions. Work in progress, tentative title. 
[27] E. Hunsicker, V. Nistor, and J. Sofo. Analysis of periodic Schrodinger operators: regularity and approximation 

of eigenfunctions. J. Math. Phys., 49(8):083501, 21, 2008. 
[28] V. Iftimie, M. Mantoiu, and R. Purice. Magnetic pseudodifferential operators. Publ. Res. Inst. Math. Sci., 

43(3):585-623, 2007. 

[29] V. Iftimie and V. Purice. Eigenfunctions decay for magnetic pseudodifferential operators, preprint. 
[30] T. Kato. Fundamental properties of Hamiltonian operators of Schrodinger type. Trans. Amer. Math. Soc, 
70:195-211, 1951. 



SCHRODINGER OPERATORS 



15 



[31] T. Kato. On the eigenfunctions of many-particle systems in quantum mechanics. Comm. Pure Appl. Math., 
10:151-177, 1957. 

[32] V. A. Kondrat'ev. Boundary value problems for elliptic equations in domains with conical or angular points. 

Transl. Moscow Math. Soc, 16:227-313, 1967. 
[33] V. Kozlov, V. Maz'ya, and J. Rossmann. Spectral problems associated with corner singularities of solutions 

to elliptic equations, volume 85 of Mathematical Surveys and Monographs. American Mathematical Society, 

Providence, PJ, 2001. 

[34] R. Lauter, B. Monthubert, and V. Nistor. Pscudodiffcrcntial analysis on continuous family groupoids. Doc. 

Math., 5:625-655 (electronic), 2000. 
[35] R. Lauter and J. Seiler. Pseudodifferential analysis on manifolds with boundary — a comparison of 6-calculus 

and cone algebra. In Approaches to singular analysis (Berlin, 1999), volume 125 of Oper. Theory Adv. Appl, 

pages 131-166. Birkhauscr, Basel, 2001. 
[36] M. Lesch. Operators of Fuchs type, conical singularities, and asymptotic methods, volume 136 of Teubner- 

Texte zur Mathematik [Teubner Texts in Mathematics]. B. G. Teubner Vcrlagsgcscllschaft mbH, Stuttgart, 

1997. 

[37] H. Li, A. Mazzucato, and V. Nistor. Analysis of the finite element method for transmission/mixed boundary 

value problems on general polygonal domains. Electron. Trans. Numer. Anal., 37:41-69, 2010. 
[38] R. Martin. Electronic structure: basic theory and practical methods. Cambridge University Press, 2004. 
[39] R.B. Melrose. The Atiyah-Patodi-Singer index theorem. Research Notes in Mathematics (Boston, Mass.). 4. 

Wellesley, MA: A. K. Peters, Ltd.. xiv, 377 p. , 1993. 
[40] R.B. Melrose. Geometric scattering theory. Stanford Lectures. Cambridge University Press, Cambridge, 1995. 
[41] S. Moroianu and V. Nistor. Index and homology of pseudodifferential operators on manifolds with boundary. 

In Perspectives in operator algebras and mathematical physics, volume 8 of Theta Ser. Adv. Math., pages 

123-148. Theta, Bucharest, 2008. 
[42] S. Moroz and R. Schmidt. Nonrclativistic inverse square potential, scale anomaly, and complex extension. 

Preprint hep-th/0909.3477v3, 2010. 
[43] C. Parenti. Operatori pseudodiffcrentiali in R n e applicazioni. Annali Mat. Pura ed App., 93:391-406, 1972. 
[44] E. Schrohe. Spectral invariance, ellipticity, and the Fredholm property for pseudodifferential operators on 

weighted Sobolev spaces. Ann. Global Anal. Geom., 10(3):237-254, 1992. 
[45] B.W. Schulze. Boundary value problems and singular pseudo-differential operators. Wiley-Intcrscience Series 

in Pure and Applied Mathematics. Chichester: John Wiley & Sons., 1998. 
[46] C. Schwab and R. Stevenson. Space-time adaptive wavelet methods for parabolic evolution problems. Math. 

Comp., 78(267): 1293-1318, 2009. 
[47] A. Vasy. Propagation of singularities in many-body scattering. Ann. Sci. Ecole Norm. Sup. (4), 34(3):313— 402, 

2001. 

[48] J. Vazquez and E. Zuazua. The Hardy inequality and the asymptotic behaviour of the heat equation with an 

inverse-square potential. J. Funct. Anal, 173(1):103-153, 2000. 
[49] H. Wu and D.W.L Sprung. Inverse-square potential and the quantum votex. Physical Review A, 49:4305-4311, 

1994. 

[50] H. Yserentant. Regularity and approximability of electronic wave functions, volume 2000 of Lecture Notes in 
Mathematics. Springer- Verlag, Berlin, 2010. 

Eugenie Hunsicker, Department of Mathematical Sciences, Loughborough University, Loughbor- 
ough, Leicestershire, LE11 3TU, UK 

E-mail address: E . HunsickerSlboro . ac .uk 

Hengguang Li, Department of Mathematics, Wayne State University, Detroit, MI 48202, USA 
E-mail address: hliamath.wayne.edu 

V. Nistor, Pennsylvania State University, Math. Dept., University Park, PA 16802, USA, and Inst. 
Math. Romanian Acad. PO BOX 1-764, 014700 Bucharest Romania 
E-mail address: nistorOmath.psu.edu 

Ville Uski, Department of Mathematical Sciences, Loughborough University, Loughborough, Le- 
icestershire, LE11 3TU, UK 

E-mail address: V.Uski@lboro.ac.uk 



